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Foldy-Wouthuysen transformation for a Dirac-Pauli dyon and the 
Thomas-Bargmann-Michel-Telegdi equation 
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^Department of Physics, Beijing Normal University, Beijing 100875, China 

The classical dynamics for a charged point particle with intrinsic spin is governed by a relativistic 
Hamiltonian for the orbital motion and by the Thomas-Bargmann-Michel-Telegdi equation for the 
precession of the spin. It is natural to ask whether the classical Hamiltonian (with both the orbital 
and spin parts) is consistent with that in the relativistic quantum theory for a spin-1/2 charged 
particle, which is described by the Dirac equation. In the low-energy limit, up to terms of the 7th 
order in 1/Eg {Eg — 2mc^ and m is the particle mass), we investigate the Foldy-Wouthuysen (FW) 
transformation of the Dirac Hamiltonian in the presence of homogeneous and static electromagnetic 
fields and show that it is indeed in agreement with the classical Hamiltonian with the gyromagnetic 
ratio being equal to 2. Through electromagnetic duality, this result can be generalized for a spin-1/2 
dyon, which has both electric and magnetic charges and thus possesses both intrinsic electric and 
magnetic dipole moments. Furthermore, the relativistic quantum theory for a spin-1/2 dyon with 
arbitrary values of the gyromagnetic and gyroelectric ratios can be described by the Dirac-Pauli 
equation, which is the Dirac equation with augmentation for the anomalous electric and anomalous 
magnetic dipole moments. The FW transformation of the Dirac-Pauli Hamiltonian is shown, up to 
the 7th order again, to be also in accord with the classical Hamiltonian. 



PACS numbers: 03.65.Pm, ll.lO.Ef, 71.70.Ej 



I. INTRODUCTION 



The relativistic quantum theory for a spin-1/2 point 
particle is described by the Dirac equation [1| . The wave- 
function used for the Dirac equation is the Dirac bispinor, 
which is composed of two Weyl spinors corresponding to 
the particle and antiparticle parts. Rigorously, the Dirac 
equation is self-consistent only in the context of quan- 
tum field theory, in which the particle- antiparticle pairs 
can be created. In the low-energy limit, if the relevant 
energy (the particle's energy interacting with electromag- 
netic fields) is much smaller than the Dirac energy gap 
Eg — 2mc^ [m is the particle mass), the probability of 
creation of particle-antiparticle pairs is negligible and the 
Dirac equation is adequate to describe the relativistic 
quantum dynamics of the spin-1/2 particle without tak- 
ing into account the field-theory interaction to the an- 
tiparticle. 

The Foldy-Wouthuysen (FW) transformation is one of 
the methods developed to investigate the low-energy limit 
of the Dirac equation ^} In the FW method, 1/Eg is 
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^ It is often said that FW method gives the nonrelativistic limit 
of the Dirac equation. The phrase "nonrelativistic" is somewhat 
misleading as it usually refers to "low-speed" limit. As we will 
show in this paper, the FW transformation (if performed to or- 
ders high enough) actually agrees with the relativistic classical 
dynamics even when the speed of the particle is large. The ap- 
propriate description is to say that the FW transformation yields 
"low-energy" limit. 



treated as the small parameter; the Dirac Hamiltonian 
in the Dirac bispinor representation is block diagonal- 
ized up to a certain order of 1/Eg and the remaining 
off-diagonal matrices, which correspond to the particle- 
antiparticle interactions, are brought into the next order 
of 1/Eg and thus neglected. This is achieved by a series 
of successive unitary transformations performed on the 
Dirac Hamiltonian. Furthermore, a series of successive 
transformation in FW method can be reduced into one 
single transformation by the use of the Lowding parti- 
tioning method Q. For a charged spin-1/2 particle sub- 
ject to a non-explicitly time-dependent field, an exact 
FW transformation has been found by Eriksen ^, and 
the validity of the transformation is studied in [5|. 

Alternatively, the Dirac Hamiltonian can also be ex- 
panded in powers of Plank constant /i [y]. In this ap- 
proach, the small parameter is not the particle's energy 
(divided by Eg), but its wave length. A diagonalization 
procedure based on the expansion in powers of h has been 
constructed in [3,[g|. In this procedure, the Berry phase 
correction can also be taken into account. Furthermore, 
the semiclassical ?i-expansion enables us to describe the 
quantum corrections on the classical expression in strong 
fields 0. 

On the other hand, the classical (non-quantum) dy- 
namics for a relativistic point particle endowed with 
charge and intrinsic spin in static and homogeneous elec- 
tromagnetic fields is well understood. The orbital mo- 
tion is govern by the relativistic Hamiltonian and the 
precession of the spin by the Thomas-Bargmann-Michel- 
Telegdi (TBMT) equation [l^. The relativistic Hamil- 
tonian for the orbital motion plus the Hamiltonian ob- 
tained from the TBMT equation (called TBMT Hamil- 



tonian) is expected to provide a low-energy description 
of the relativistic quantum theory. The conjecture that 
the low-energy limit of the Dirac Hamiltonian reduces to 
the classical orbital Hamiltonian plus the TBMT Hamil- 
tonian has been suggested but remains to be affirmed. 

In order to investigate the consistency between the low- 
energy limit of the Dirac equation and the classical dy- 
namics, we perform a series of FW transformations and 
expand the Dirac Hamiltonian up to terms of the 7th 
order in 1/Eg. The electromagnetic fields are assumed 
to be static and homogeneous. Taking care of the rela- 
tion between the kinematic momentum used in the Dirac 
Hamiltonian and the boost velocity used in the TBMT 
Hamiltonian, we show that the FW transformation of 
the Dirac Hamiltonian is in agreement with the classi- 
cal orbital Hamiltonian plus the TBMT Hamiltonian for 
the case of the gyromagnetic ratio equal to 2. Through 
electromagnetic duality, this result can be generalized for 
a spin-1/2 dyon [ll|, which has both electric and mag- 
netic charges and thus possesses both intrinsic electric 
and magnetic dipole moments (with both gyromagnetic 
and gyroelectric ratios equal to 2). 

To affirm the consistency to a broader extent, we need 
to show that the relativistic quantum theory of a spin- 
1/2 dyon with arbitrary values of the gyromagnetic and 
gyroelectric ratios also reduces to the classical counter- 
parts as a low-energy limit. The relativistic quantum the- 
ory of a spin-1/2 dyon with the inclusion of anomalous 
magnetic dipole moment (AMM) and anomalous electric 
dipole moment (AEM) can be described by the Dirac- 
Pauli equation [9|, |l2[ , which is the Dirac equation with 
augmentation for AMM and AEM. The FW transforma- 
tion is performed on the Dirac-Pauli Hamiltonian, again 
up the 7th order in 1/Eg, and the result confirms that it 
remains in agreement with the classical orbital Hamilto- 
nian plus the TBMT Hamiltonian for arbitrary values of 
the gyromagnetic and gyroelectric ratios. 

This paper is organized as follows. In Sec. |TT1 we in- 
vestigate the tensorial structure of the orbital and intrin- 
sic dipole moments. In Sec. IIIIl we briefiy review the 
classical orbital Hamiltonian and the TBMT equation. 
In Sec. IIVI we perform the FW transformation on the 
Dirac Hamiltonian for a spin-1/2 dyon and show that it 
agrees with the TBMT equation for the case of the gy- 
romagnetic and gyroelectric ratios equal to 2. Later in 
Sec.|Vl we perform the FW transformation on the Dirac- 
Pauli Hamiltonian and show that it again agrees with 
the TBMT equation even with the inclusion of AMM 
and AEM. Finally, the conclusions are summarized and 
discussed in Sec. |VT1 Some calculational details are sup- 
plemented in Appendices [K\ and IB] 



II. ORBITAL AND INTRINSIC DIPOLE 
MOMENTS 

For a general Lorentz transformation from the primed 
(boosted) frame to the unprimed (laboratory) frame, the 



transformation of a 4- vector k'^ = (k" , k) is given by J13I : 

(1) 



k0 = 7(k'0+/3-k'). 



k = k' 



/32 



(/3-k')/3 + 7k'"/3, 



where v = c/3 is the boost velocity of the primed frame 
relative to the unprimed frame, 7 is the Lorentz factor 
7 = l/yT^Fand/3=|/3|. 

In the primed system, let us consider the case that 
charge and current densities satisfy the conditions: 



d'xp = 0, 



(fx'3' = 0. 



(2) 



The vanishing of the total charge means that the system 
is neutral^ and the vanishing of the total current is a 
consequence of the static condition: dp' /dt' = — V'- J' — 
Q? Because the charge density and current density form 
a 4-vector J^^ = {cp,3), it can be shown that the same 
conditions also hold in the unprimed system: 



V Jv 



(3) 



In the unprimed frame, the magnetic dipole moments 
m is defined as 



where 



m = / d^Xflrn, 

Jv 



Zc 



(4) 



(5) 



is the magnetic dipole moment density, and the electric 
dipole moment p is defined as 



where 



p = 2 / d^xti'p, 
Jv 



Mp = 2^P 



(6) 



(7) 



is the canonical electric dipole moment density (the extra 
factor of 2 is introduced for later convenience). In the 
primed system, the definitions of both dipole moments 
are the same as those in the unprimed system. Using 
Eq. (H)), Eq. ^ can be written as 



7 - 1 

Mm =Mm + -Ip-f^ ^ (Mm >< /3) 
+ i7(x'cp' - X'^J') X /3, 



(8) 



^ Since the static condition gives V' • J' = 0, it can be shown J^ ■ 
V ■ (x^J'). Consequently, fy,d'^x'J[ = J^,(fx'SJ' ■ {x'^3') ■■ 
Jgyi da! ■ {x'^3') = if the current J' is locaUzed. 



where /i.'„ = x' x J'/2 is the magnetic dipole density in 
the primed system. It is interesting to note that if we 
integrate the term x'cp' — a;'° J' in the primed system, we 
obtain 



/ ax [x cp — X J ) = / ax x cp = cp , 
Jv Jv 



(9) 



where the neutral condition [Eq. ^] has been used. This 
suggests that we can define the tensorial electric dipole 
moment as 



fJ'P 



2c ^ 



a;"J 



(10) 



so that the dipole moment can be defined as a second 
rank antisymmetric tensor: 



M'"' 



1 
2^ 



{x^r -x^J^"). 



(11) 



The canonical and tensorial dipole moments densities 
yield the same (integrated) dipole moments, because the 
neutral condition ensures that the integration of the sec- 
ond term x^3 vanishes in the unprimed system. The 
components of the second rank tensor iW"^ are 



M" 



l-y^^' 



x\f) 



M'^ = —{x\P-xU') 



-Mp, 



eijkp.„ 



(12) 



Consequently, the Lorentz transformation between 
(/Xp, Urn) and {fip, ^i'„^) is of the form 



Mp = 7(Mp + /3 X /^m) rT/3(/3 • Mp), 

7 + 1 



(13) 



Mm = 7(Mm - /3 X Mp) 



7+1 



/3(/3-Mm)- 



The transformation [Eq. P^ ] is exactly the same as that 
for the electric and magnetic fields if we take the replace- 
ment rules: E o /Xp Bo —fim- 

The corresponding transformation for the (integrated) 
dipole moments is given by 



P 

2 ' 

m 



— + /3 X m' - 

2 ^ 7+1 



' m-\) 



7 



m 



-/3x^-^- 

2 7+1 



/3(/3-m') 



(14) 



where dPx = j(Px' is used. Since the extra factor 7 
arises in the right hand side of Eq. ([T^ due to the spatial 
integral, p/2 and m do not transform covariantly and 
thus do not form a second rank tensor unlike fip and 

In the case with only proper electric dipole moment 
and no proper magnetic dipole moment (i.e. p' 7^ and 
m' = 0), when boosted, the electric dipole will result in 
a magnetic dipole moment m = — 7^(/3 x p'/2) in the un- 
primed system. This can be understood as follows: If we 



think p' as two endpoints separated by a short distance 
and charged with +q and —q, in the unprimed system, 
the positive and negative charges acquire a velocity and 
give rise to currents in opposite directions, thus result- 
ing in a magnetic dipole moment. In an inhomogeneous 
magnetic field, a moving object with only proper electric 
dipole moment can feel the magnetic force F = (m- V)B. 

On the other hand, in the case with only proper mag- 
netic dipole moment and no proper electric dipole mo- 
ment (i.e. m' 7^ and p' — 0), when boosted, the 
magnetic dipole will result in an electric dipole moment 
p = 27^ (/3 x m') in the unprimed system. This is due 
to the fact that in the unprimed system the charge den- 
sity p arises through the Lorentz transformation even if 
the charge density is zero in the primed system {p' = 0). 
The charge density in the primed system cp = 7/3 • J' 
is positive (negative) when the current is parallel (anti- 
parallel) to the boost velocity; therefore, as a magnetic 
dipole can be thought as a small current loop, the small 
current loop in the primed system gives rise to opposite 
charges separated by a short distance in the unprimed 
system, thus resulting in an electric dipole moment. In 
an inhomogeneous electric field, a moving object with 
only proper magnetic dipole moment can feel the electric 
force F = (p • V)E. 

The dipole moments considered above are orbital in 
the sense that they are sourced by the orbital distribu- 
tion of J^(a;). On the other hand, a point particle can 
give rise to an intrinsic dipole moment if it is charged 
and endowed with intrinsic spin. The fact that fip and 
Hm form an antisymmetric tensor iW"^ suggests that the 
intrinsic spin s can be generalized to a second-rank an- 
tisymmetric tensor S^'^, which gives the intrinsic dipole 
moments as 



2mc 



(15) 



where e is the electric charge of the particle, m the mass 
and ge the gyromagnetic ratio. The spin has only three 
independent components; thus S'^'^ is dual to an axial 
4- vector 5" = (5°,S) via 



and conversely 



S"" = - e''''°'^U^Sfi 



^" = ^ e^^^%S^s, 



(16) 



(17) 



where U°' is the particle's 4-velocity. The 4-vector S" 
reduces to the spin s in the particle's rest frame; i.e., 
gfa ^ (5"0^s') = (0,s). The vanishing of the time- 
component in the particle's rest frame is imposed by the 
covariant constraint: 



(18) 



C/„S'" = 0. 



In the particle's rest frame, C/'" — (c, 0, 0, 0) and Eq. p? 
yields 



/ _ gee 
'^" " 2mc ^' 



m; = o. 



(19) 



Therefore, the intrinsic spin gives only the proper intrin- 
sic magnetic dipole and no proper intrinsic electric dipole. 
In order to have both proper intrinsic magnetic and elec- 
tric dipoles, we consider a dyon particle [ll!|, which pos- 
sesses both electric charge e and magnetic charge e, and 
Eq. (fTSJ) is generalized as 



M"" = M'"' + Air = -^ S"" + ^^ S"", (20) 
2mc 2mc 

where gg is the gyroelectirc ratio and 






(21) 



is the dual of S^"" . In the rest frame, Eq. (gOl) yields both 
magnetic and electric dipoles: 






2mc 



2mc 



(22) 



III. THE 

THOMAS-BARGMANN-MICHEL-TELEGDI 

EQUATION 

Consider a relativistic point particle endowed with 
electric charge and intrinsic spin subject to static and 
homogeneous electromagnetic fields. The orbital motion 
of the particle is described by 



"dT 



mc 



(23) 



and the precession of the spin is govern by the TBMT 
equation [10|: 



"dT 



e 
mc 



■ F^f^SB 



1 



9e 



2- -''5 + ?lY-i)^"('^^^'"^A<) 



(24) 

Equation (|23p in the covariant form can be shown to 
be equivalent to the Hamilton's equations : 



dx 
dp 



— {P,-fforbit} 



(25) 



in the unprimed frame, where p is the conjugate mo- 
mentum to X and the Hamiltonian i?orbit governing the 
orbital motion is given by 



iTorbit (x, p) = v^(cp-eA(x))'+m2c4 + e 0(x) (26) 

with A" = {(/), A) being the 4-vector potential for the 
electromagnetic field F^^"^ . (See Sec. 12.1 in fl3'| for more 
details.) 

On the other hand, Eq. ([M)) leads to 



ds 



mc 



■ s X F(x) 



(27) 



with 



P.,|_:,iB-(|-i 



7 
7 



7+1 



(/3 • B)/3 



7+1 



(28) 



/3xE, 



which gives the spin precession with respect to the time 
of the unprimed frame. (See Sec. 11.11 in |l3| for more 
details.) Because {si, Sj} = eijkSk, Eq. (|27p can be recast 
as the Hamilton's equation: 



ds 
di 



\S, J^spinj 



with 



H. 



spin(x,s) = S-F(x) 

mc 



(29) 



(30) 



called the TBMT Hamiltonian, which governs the pre- 
cession of the electric dipole subject to a static and ho- 
mogeneous field. 

In the low-speed limit {j3 ^ 1), we have 7 w 1 and 
Eq. (1201) gives 



ffs 



2mc 



5eB - (I - 1) (/3 . B)/3 



-(ffe - l)/3 X E 



(31) 



The first term in Eq. (|31|) is the interaction energy of 
the magnetic moment /x^^ in the magnetic field, which 
accounts for the anomalous Zeeman effect. The second 
term corresponds to the change rate of the longitudinal 
polarization, which vanishes in the case of g^ = 2. The 
third term is the spin-orbit interaction (the interaction 
of the boosted electric dipole fip ~ f3 x fi'^ coupled to 
the electric field) plus the correction for the Thomas pre- 
cession. 

By treating x, p and s as independent phase space 
variables, the total Hamiltonian is given by'^ 

H{X, p, S) = i7orbit(x, p) + -ffspin(x, s). (32) 

If the particle has both electric charge e and magnetic 
charge e (i.e. the particle is a dyon), Eq. ^^\ and 
Eq. (|30p are modified with the inclusion of the dual coun- 
terparts; i.e. 



-fforbit(x, p) = -\/ fcp-eA(x) -eA(x)j +m'^c^ 

' e(/i(x) + g0(x) (33) 



^ Note that, in order to add -ff orbit and Hgpin together, we have to 
consider ds/dt = dS' /dt in Eq. II27I I. instead of dS/dt, dS/dr or 
ds/dr. This is because Si are degrees of freedom independent of 
X and p, but Si are not. Furthermore, to be consistent with the 
orbital motion, the precession is cast with respect to t, instead 
of the proper time r of the moving particle. 



and 



F,pi„(x,s) = -— s-F(x)-— s-F(x) (34) 

mc mc 



with 



■'=lf-l + ^)B-(f-l)^(/3^W 



5e 7 



2 7 + 1 



(35) 



/3xE, 



where A — {(f), A) is the dual 4- vector potential which 
gives Pf^" = df'A'' - d^'A'^ and F'^'^ := l/2e'^''"^Fa/3 is 
the dual field strength (i.e. B = -E and E = B). 

Equation (1^ and the TBMT equation given in 
Eq. (P^ are derived as the requirement of covariant is 
considered. They are classical (non-quantum) equations 
and we wonder whether the Hamiltonian given in Eq. (j32p 
is consistent with that in the relativistic quantum theory 
for a charged point particle with intrinsic spin. 

The relativistic quantum theory of a spin- 1/2 particle 
is described by the Dirac equation. The Dirac bispinor 
however has both the particle and antiparticle compo- 
nents, which are entangled by the Dirac equation. 

In order to compare with the TBMT equation, we con- 
sider the low-energy limit in which the relevant energy 
is much smaller than the Dirac energy gap Eg and the 
FW transformation is used to block-diagonalize the Dirac 
Hamiltonian. In Sec. IIV[ we will show that the FW 
transformation of the Dirac Hamiltonian indeed agrees 
perfectly with the TBMT equation up to the 7th order 
of 1/Eg with the intrinsic spin given by s = ha/2 [ui 
are the Pauli matrices) and the gyromagnetic ratio given 
by ffe — 2. This can be easily generalized for a Dirac 
dyon by adding the magnetic charge (and we will have 
5e := gg = 2). 

As the Dirac equation always yields ^e = 2, we will 
not see the second term in Eq. ((28|). which accounts for 
change of the longitudinal polarization. In order to see 
that the quantum theory is in accord with the TBMT 
equation even for the case of ^e y^ 2, we study the Dirac- 
Pauli equation in Sec. |V]with the inclusion of anomalous 
dipole moments. The results again affirms the consis- 
tency between the FW transformation of the Dirac-Pauli 
Hamiltonian and the TBMT equation up to the 7th order 

of l/Eg. 



IV. FOLDY-WOUTHUYSEN 

TRANSFORMATION FOR THE DIRAC 

HAMILTONIAN 

The relativistic quantum theory of a Dirac particle is 
described by the Dirac equation 



ih^^\^)^H\i,), 



(36) 



where the Dirac bispinor |'0) = (x, Lp)^ is composed of 
two 2-component Weyl spinors x ^-nd (/? corresponding to 



the particle and antiparticle parts. The Dirac Hamilto- 
nian is given by 



H = mc^P -I- c a • n + y^, 



(37) 



where the 4x4 matrices /3 = CTz (8) 1 and ai = ax® en are 
given in the Pauli-Dirac representation ll| and satisfy** 



{/3,d.} = 0, 
{di,dj} = 2(5j. 
d? = /32 = 1. 



(38) 



The mass of the particle is m and 11 is the kinetic mo- 
mentum. 

In order to have the proper intrinsic electric dipole 
moment and the proper intrinsic magnetic dipole moment 
at the same time, we consider a Dirac dyon (i.e. a Dirac 
particle with both electric charge e and magnetic charge 
e). For a dyon }lli] . the kinetic momentum is given by 



n = p -A -A, 

c c 



(39) 



and the scalar potential V is composed of electric and 
magnetic monopole potentials: 



V = e4> + eq 



(40) 



In the following, we will first perform the successive 
FW transformations of the Dirac Hamiltonian up to the 
7th order of l/Eg in Sec. IIVAI and later show that the 
results agree with the TBMT equation in Sec. IIVBI 



A. Foldy-Wouthuysen transformation 

In order to perform the FW transformation [3] , we have 
to rewrite the Dirac Hamiltonian to the form: 



2 



H = Eg^ + no + nE, 



(41) 



where Eg = 2mc^ is the Dirac energy gap, and the odd 
matrix flo and the even matrix ^Ie are defined as 

{/3,f}4 = 0, [/3,J7b]=0. (42) 

In the case of Eq. ([37]), we have 

no^ca-U, nE = V. (43) 

The resulting effective hamiltonian iJpw can be ob- 
tained by the successive unitary transformations which 
partitioning off the odd matrices to a higher order. In 



^ To avoid confusion with the boost velocity /3, we use the checked 
notation /3 to denote the 4x4 matrix. With the same style, the 
notations a, ■y [defined in Eq. I|88| l] and fj [defined in Eq. (|95| l] 
are checked as well. 



general, we can use the FW matrix t/pw as a single trans- 
formation, and expand the exponent of the matrix in 
powers of 1/^; this is the well-known Lowdin partition- 
ing method li\. It can be shown that the FW transfor- 
mation of Eq. (|4ip . namely the transformed Hamiltonian 
denoted as iJpw — UywHU^^, up to terms of the 7th 
order in 1/Eg can be written as (see Appendix |A|) 



H 



FW 



PEg 



^E + Y.H'il. + oil/El), (44) 



1=1 



M 



where the first four terms H^y^, I = 1,2,3,4 are given 
by 



B, 



(1) _ ml 



FW 



^. ' 



„(2) _ ^ / W 

™ ~ El\2 

rrO) _ 
^FW - p3 
9 



(45a) 

(45b) 

Li^^pnl + P{pvf], (45c) 

<w = ;^ (^[[^o, W]f}o] - \[DMl]) , (45d) 

and the operators W and V are defined as 

P=[fio,f}£], (46a) 

W=[V,Vlo]. (46b) 

By using Eq. (|43p . it can be shown the three Hamil- 
tonians -ff^-w ' ' are in agreement with the previous re- 



sults 



a a 



w 
. The term of the 5th order is given by 

3 



e,m=i 

(£+m=4) 



(£+m+n=3) 



(47) 



where the subscript (5) in the commutator [{$fto)(5), ^o] 
indicates that the commutation of /351o with flo is per- 
formed successively by five times; i.e., [{$no){5),^o] = 

[prio, [(3^o, [(3^o, [$no[(3no, ^o]]]]]. The term of the 6th 
order is 

e,m=i 

(e+m=5) 

+ 1 E E[/30''M/3o<"\/^<"^]], 



e,m=ln=Q 
(l+m+n=A) 



where the odd matrices O*^^^ for l — 1, 2, 3, 4 are given by 



6 15 

and the even matrices /i*^"-' for n = 0, 1, 2 are 



(49) 









/.(I) 



P^l 



(50) 



rW 



To obtain the FW transformed Hamiltonian H^^ up to 
the 7th order of I /Eg, we need only three successive 
transformations C/pw = exp(S'3) exp(S'2) exp(S'i) (see 
Appendix IX)) . and it can be shown that Si, S2 and S3 
are all anti-hermitian matrices. 

To simplify the calculation, some restrictions and as- 
sumptions are made. The electromagnetic field is as- 
sumed to be static, as has been used in obtaining 
Eq. (|1^ . In order to demonstrate the equivalence clearly 
between the TBMT Hamiltonian and -ffpwj we further 
assume that the external fields E and B are homoge- 
neous, and thus the field gradient vanishes. Furthermore, 
the terms proportional to products of field strengths, 
such as EiEj, EiBj and BiBj, are all neglected as a 
good approximation for weak fields. 

We now evaluate each term of B^t^^. The kinetic term 
^1/Eg in Eq. (|45ap can be written as 



fjg (eg ■ n)^ 

Eg 2?77,C^ 

2m 



(51) 



where [on, cnj] — 2ieijk'Sk is used. By using the definition 
of magnetic field B = V x A and dual magnetic field 
B = V X A, we have cII x 11 = ih{eB + eB). By 
applying the duality B = -E in Eq. ^T^, Eq. P5a|) then 
gives 



B, 



(1) 

FW 



mi 

2m 



-/5 



eh 
2mc 



S .B~/3 



eh 
2mc 



S -E. (52) 



The second term of Eq. ([5^ is the Zeeman Hamiltonian 
for an electron (with e = — |e|) [TSJ, and the third term is 
its duality. It is interesting to note that e/iS/2mc plays 
the role of electric dipole moment because it couples to 
the electric field. In this sense, we can define the (proper) 
intrinsic electric dipole moment (n'p) and (proper) intrin- 
sic magnetic dipole moment (/x^) as 



M« 



Mp 



— S 
2toc 

eh ^ 



2mc 



(53a) 
(53b) 



Equation (j53al) implies that the dyon's intrinsic gyromag- 
(48) netic ratio ge = 2 for the Dirac Hamiltonian. We also find 



2 for the dyon's intrinsic 



r(2) 



the same gyroelectric ratio g^ 
electric dipole moment. 

We now focus on the 2nd-order Hamiltonian H'^^^ m 
which the spin-orbit coupled term is included. It can be 
shown th at W is given by W = -2c^hi: ■ ((eE + eE) x n) 
and Eq. (|45b[) can be written as 



H, 



(2) 
FW 



IE 
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X (J'r. 



1. 



-B 
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X/^p 



(54) 



where the duality E = B is used. The first term of 
Eq. (|54l) is the spin-orbit interaction for an electron |15| . 
On the other hand, we neglect the terms proportional to 
products of field strengths in evaluating the two terms in 
Eq. (|45cl) . and thus we can obtain 



H, 



(3) 
FW 






(M™ • B) 



(55) 



E) 



where the assumption of homogeneous electromagnetic 
fields is used, and thus the operator |np commutes with 
the magnetic field B. If the magnetic charge e vanishes, 
the first term of Eq. ([55]) is the relativistic mass correction 
that contributes to the spectrum of fine structure 15]. 

The second term of Eq. (|55|) is the relativistic cor- 
rection to the Zeeman Hamiltonian appearing in Hp^ 
[Eq. ((52)) ]. For the 4th order H^tJ^, it can be shown that 



24 



[K,>V],f^o 
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(56) 
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by a relativistic correction — 3(|n|/mc)^/4. For those 
terms in the 5th order, it can be sown that each term 
corresponding to Eq. (|47| is given by 

"* 1 QO 

^ [/30W,0('")] = --/3{2?,[f7„>V]} + i-/3176, 

2 1 

{e+m+n=3) 

7 18 

= ~^P{v, [no, w]} + —{V, Qovno}. 

(59) 



We note that the operator V is proportional to a • E, 
which is of the 1st order of the electric field as well as the 
operator W. We find that the terms {V, ['[lo,W]} and 
{VjQoT^no} in Eq. (|59l) are proportional to the product 
of only electric field, and thus will be neglected. On the 
other hand, the magnetic field in f]^ is also of the 1st 
order [see Eq. (ISTj) ]. If we further neglect those terms 
proportional to the products of magnetic field and con- 
sider the homogeneous field, ft^ becomes 
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If we neglect all terms proportional to products of elec- 
tromagnetic fields, one can show that (sec Appendix IBJ) : 



Therefore, H^.^. with substitution of Eqs. §^ and (pO)) 
becomes 



noW^o 



2c^\n\^w, 



c^inj^w. 



(57a) 
(57b) 



Note that there is a minus sign in Eq. (j57b[) . The 4th- 
order term Eq. (j45dp with substitution of Eqs. (|56| and 
(1571) can be written as 
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It is interesting to note that the 4th order Hamiltonian 
TJp^ is in relation to the 2nd order Hamiltonian H^^ 



The first and second terms of Eq. (|6T|) contribute to the 
relativistic mass correction and the Zeeman effect, re- 
spectively. For the 6th-order term Hp^, it can be shown 
that the commutators of the form [/^O^, O*^™^] are given 
by 






(62) 



where we also have [/30(3), 0(2)] :^ [/30(2), ©(3)] and [/30(4),0(i)] = [l30^^\0^^^. On the other hand, the commutators 



of the form [/SO^^), [/SO^"), /i^")]] in iJ*,''^ are given by 



Ifi 
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[^0(1)^ [^0(2)^ /,(!)]] ^ ^(r)4H; + r)3n;r!o + nlwnl + noWnl + wnt), 



(63) 



Because W is of the 1st order of an electric field as well 
as P, we can use Eq. (I57p to reduce these equations into 
a form with only fields of the 1st order. For example, 
the term nlwrio becomes nlVVflo = nHQoWno) 



teraction in Hj^^. Therefore, H^^ and Hp^ and -ffp^y 
are composed of kinetic energy, interaction energy of Zee- 
man effect and their relativistic corrections. On the other 
hand, Hp^ and Hp^ and i^p\y contain only spin-orbit 



|np(— |npyV). On the other hand, [P, [P,W]] and interaction and its relativistic corrections. 



[D, [[flo, W], ^e]] are neglected because they are at least 
of the 2nd order of fields. In that sense, by the use of 
Eqs. dnH and (p5)) . one can obtain 



To simplify the expression of ffpi^, we can define a 
scaled kinetic momentum operator ^ as^ 



[(/317o)(6),f]s]«16c4|n|4w, 
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J2 Yl [/3C^^\ [/3C^™\ ^^"^]] ~ ^c^ini'^w. 



(£+m+n=4) 
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Therefore, Eq. P5|) with substitution of Eq. (IMl) becomes 
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Eq. (1651) is the relativistic correction to the spin-orbit in- 



n 

mc 



(66) 



By replacing |n|/r7ic with Eq. (|B5)) . Eq. dill) with sub- 
stitution of Eqs. ([521) . ([541). ([FS ]) . ([55 | . (|6T1) and (|65l) be- 
comes a sum of two terms: 



7?FW ~ ^orbit + Hgpin, 



(67) 



where -fforbit is the kinetic energy plus the potential en- 
ergy, namely 



iJorbit - ^mc^ ( 1 + lie - l^^ + ^m' ] + V, (68) 



and i^spin is the energy of intrinsic dipole moments plac- 
ing in electromagnetic fields, namely. 



^ It must be stressed that the operator ^ does not directly cor- 
respond to the Lorentz boost velocity /3 given in the previous 
section. The appropriate transformation between ^ and opera- 



tor for the boost velocity is considered in Sec lIV Ej] 
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In Sec. IIVB| we will focus on the dipole Hamiltonian 
[Eq. (|M|) ]. We will show that Eq. ([S^ is in agreement 
with TBMT equation, provided that the proper transfor- 
mation of ^ and Lorentz boost velocity (3 is taken care 
of. 



B. In relation to TBMT equation 

We will show that the FW transformation of the Dirac 
Hamiltonian of a dyon is equivalent to the Hamiltonian 
obtained from TBMT equation with g^ = g^ — 2. That 
is, Eq. dSH) is equivalent to Eq. ^ and Eq. ^ to 
Eq. dill) with Eq. ^ and Eq. ^ for g^ ^ g^ ^ 2. 
However, we must first find the boost velocity in order to 
compare them. It must be emphasized that /3 in TBMT 
equation is the boost velocity but ^ in -ffpw is not. One 
has to define the boost operator (3 via 



because the kinetic momentum 11 = mc$, = m\J and 
the 4- velocity U"' = (70,7^). By using Eq. ([70]) . 
the kinetic energy operator Eq. (|68| behaves like 

2(. , 1,^|2_1|^|4 + ^|^|4) ^ ^,2(1+ 1|3|2 + 

0(8)). On the other hand, the expansion 
13"^ with respect to small 



mc" (1 

of Lorentz factor 7 



H(i?- 



l/v^ 



boost /3 is 7 = 1 -h i/32 + f/?-* -h fg/?** + 0(8). This im- 
plies that the the kinetic energy operator corresponds to 
the classical relativistic energy 7mc^, as expected. The 
boost operator /3 plays an important role on showing the 
equivalence between iJspin and TBMT Hamiltonian. For 
an electron, Eq. ((M)) without a magnetic charge (e = 0) 
becomes 
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The first term in the right hand side of the first equality 
of Eq. (|7ip is an effective electric dipole moment caused 
by the boosted intrinsic spin magnetic moment, which is 
the spin-orbit interaction. The second one is the Zeeman 
term. Nevertheless, Eq. (|71l) provides the relativistic cor- 
rection to the Zeeman and spin-orbit interactions. For 
the Zeeman term, the non-relativistic limit up to 1/mc 
is 



^, the relativistic correction to -ffzccman is 



Fznnn... = -/3 ( 1 " ^l^ + ^1^ ] M™ ' B. (73) 



On the other hand, the spin-orbit interaction denoted as 
Hso is 



-(^Kn ■ B, 



(72) 



which is the same as the interaction of a classical mag- 
netic moment and a magnetic field. To the 4th order of 



H„, 



-An • I X E 
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where e = — |e| is used in the second equahty, and the 
relativistic correction to this term is 



Hff, 



^" 









Kn ■ I X E. 



(75) 



We now go back to the discussion of H^^^^ and TBMT 
equation. In order to compare Eq. (fTTj) with TBMT 
Hamiltonian, we have to transform ^ in Eq. ([71]) to /3. 
Using Eq. ([70]) . we have 



2 IS 



k\' 



i_a!_ffi+„,e). (76.) 



'^-|«p 



111 = 



m i/3i= 



16 



-0(7). 
(76b) 



Therefore, we show that up to the fifth order of boost 
velocity /], the Dirac Hamiltonian of an electron is equiv- 
alent to the TBMT Hamiltonian which is obtained from 
the requirement of covariance form of classical spin. This 
implies that in the FW representation, after summing 
over all infinite expansion terms, the Dirac Hamiltonian 
of an electron would be of the form 



h: 



e=0 



fJ-n 



filB 

, 7 



1+7 



for the spin part, and of the form 



^orbit = ^Pmc + V 



/3 X E 



(80) 



(81) 



The effective spin magnetic moment in TBMT Hamilto- 



nian [Eq. ([50)1 ] transforms like {l/j)fj,'^, and we have 
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2 



P' 



0(6), 



(77) 



which is exactly the same as Eq. (|76ap up to terms of the 
4th order in /?. On the other hand, the effective electric 
dipole moment transforms like (^e/S ^7/(7+!)), and ge 
factor in the Dirac Hamiltonian is always 2. We obtain 



for orbital part. The effective magnetic field in Eq. ([50)1 
is the same as Eq. ([30]) with Eq. ([28]) for ffe = 2. 

Furthermore, for the FW transformation of the Dirac 
Hamiltonian [Eq. (l67l) ]. the TBMT equation can be gen- 
eralized to include an effective spin magnetic moment 
resulting from the boosted intrinsic electric dipole mo- 
ment. To the 1st order in j^j = |n|/TOC, we find that the 
effective dipole moments transform like 
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1 + 7 
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P P' 



16 



+ 0(7), (78) 



which is exactly the same as Eq. (I76b[) up to terms of the 
5th order in /3. Because ge factor equals 2, the longitudi- 
nal term /xj^ ■ ^ disappears in both TBMT Hamiltonian 



and H, 



(e=0) 



spin 



In the following, we will use the following 



two approximations directly; 
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(/xp,ff « p^l'; + - (I X /x;^j , 



(82) 



This means that an intrinsic electric dipole moment can 
result in an effective magnetic dipole moment when it is 
moving. Nevertheless, a moving spin magnetic moment 
can also intrinsically induce an effective electric dipole 
moment. Consider higher orders of the boost velocity, 
we rewrite Eq. (j69p as 
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(83) 



It is shown that intrinsic dipole moments transform like 
(1 - \\i\^ + 1 1^1^) w 1/7 and the boosted dipole mo- 
ments transform as 5 (1 - j\i? + fill"*) « (1-^) (see 
Eq. ([75]) ). Therefore, (/x^ )cfF and (/.tm)cff do not form a 
second rank tensor in the sense that their transformation 
is not a covariant form like Eq. (IT5)) . but the following 



form: 



(/^^)eff « pilJ^l - f 1 - ^ 

7 V 7+1 



/3 X /x: 



le 
m I ' 



f3x fi'. 



(84) 
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This implies that an energy caused by dipole moments 
in the description of Dirac Hamiltonian is not simply the 
contraction of tensorial dipole density and field tensor: 
ffspin 7^ ~fJ'p • E — fim • B, in which /Xp and /i„i transform 
as in Eq. (|T3l) . As a result, if spin is not a Lorentz scalar. 
In short, in this section we have shown that up to terms 
of the 7th order in 1/ Eg, the FW transformation of the 
Dirac Hamiltonian of an electron is in agreement with 
TBMT Hamihonian [Eq. ^] with 5e = 2. The result 
can be generalized to a particle with an intrinsic electric 
dipole moment. Because of the duality of electromag- 
netic fields, a Dirac dyon would manifest this feature. 
Furthermore, we also find the relativistic corrections to 
the Zeeman term and spin-orbit interaction 



V. FOLDY-WOUTHUYSEN 

TRANSFORMATION FOR THE DIRAC-PAULI 

HAMILTONIAN 



In Sec. IIVI we have shown that, up to the 7th order in 
I /Eg, the FW transformation of the Dirac Hamiltonian 
for a dyon is in agreement with Eq. (pS)) and Eq. ([27]) 
for 5e = ffe = 2. Since the Dirac Hamiltonian auto- 
matically yields ge — 2 and g^ = 2, the second term 
in Eq. (pS)) and Eq. psp vanishes and thus the longi- 
tudinal polarization does not change. In order to see 
that the relativistic quantum theory of a spin- 1/2 parti- 
cle is in accord with the TBMT equation even when the 
change rate of the longitudinal polarization is concerned, 
we have to study the spin-1/2 particle with anomalous 
magnetic dipole moment (AMM) and anomalous electric 
dipole moment (AEM). 

The relativistic quantum theory of a spin-/2 dyon with 
the inclusion of AMM and AEM can be described by the 
Dirac-Pauli equation [3, [l3 | 



ih^^\^) = n\^), 



(85) 



where the Dirac-Pauli Hamiltonian H is the Dirac Hamil- 
tonian H [given in Eq. ([37]) ] augmented with the correc- 
tions for the AMM and AEM: 

H = i7 + /^'(-/3S-B + i7-E) + d'(/3S-E + i7-B). (86) 

The coefficients /i' and d! are defined as follows 

3^ ,, / Qp \ eh 
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9e A ell , fgg \ 
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(87) 



which measures the AMM and AEM, respectively (note 
^' = for Qe = 2 and d' = for g^ = 2). The 4x4 
matrices /3S and 7 are defined as 



/3S 



cr 

-cr ; ' 
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cr 
-cr / ' 



where a = {(Jxi'^yT'^z) a-re Pauli matrices. We will see 
that the Dirac-Pauli Hamiltonian given in Eq. ((55)) is 



compatible to generic values of g^ and g^ and thus can 
accommodate AMM and AEM. 

In order to obtain the FW transformation of Eq. (|86p , 
we have to rewrite Eq. (|86p in terms of odd and even 
matrices. According to Eq. (|42]) . we have 



n = Pmc^ 



^E + ^t 



(89) 



where the superscript A indicates the inclusion of anoma- 
lous dipole moments. We note that /i' and d' are of the 
1st order of 1/Eg. Therefore, fJ^ and Q^ can be written 
as 



(90) 
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= Egfl', d" = Egd'. 



(91) 



The superscript / indicates that these terms are of the 
1st order of electromagnetic fields. Because we consider 
only those terms proportional to the 1st order of fields, 
the products of fields will be neglected. The validity of 
Eq. (jA17p is still true provided that the odd term of the 
second FW transformation denoted as C in Eq. ([M)) 
starts from ^/Eg. This can be seen as follows. After 
the first FW transformation. Si = (3^1^ /Eg, Eq. 
becomes 
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E„ + h + O, 
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where h and O are even and odd terms, respectively. It 
can be shown that the odd term O can be written as 
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where the corresponding 0(") , n = 1, 2, 3, 4, are given by 

o'-'^ = ^[no,nE], 

0(2) = 2c^fy(-^"B + d"E) ■ n - ^nl, 

4 1 (94) 

o(3) = __ (Qoin^, nfj + n{,nl) + -pin^, w], 
0(4) = _4 3^^(_^.B ^ ^.E) . n|n|2 + ±ni 
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where the matrix f) is defined as 

-1 

1 



(95) 



If ge and gg are both equal to 2, then Eq. (|94p goes back 
to Eq. (|49p . Using the second FW transformation S2 = 



12 



pO/Eg on Eq. ([93]) . we can obtain the other odd term 
denoted as O' . The 1st order term O'^^^ is zero because 
O starts form the 1st order of 1/Eg at least. O'^^) can 
be written as C^^^ = J/SO^^', /i(°)]. Nevertheless, we have 
0(1) = p[no,nE\ = piha. ■ (eE + gB) and /i(°) = Q-e = 
y is a scalar that commutes with 0'''^\ and thus C'(^) 
vanishes. This implies that Eq. (jA17p is still valid in this 
case. The matrix h in Eq. ([M|) is given by 
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h=ni + ^[^ntn^] + ^m^)i2),ni] 



^[(/3f7^)(3),f7^] + ^Pfi^)(4),fi^] 
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To obtain h' , we need extra corrections to h, as shown in 
Eq. (jAlSp . The resulting Hamiltonian can be written as 



T^FW — ^FW + ^FWl + ^FW2- 



(97) 



The Hamiltonian i^FW is given in Eq. (|67p . 'Hfwi con- 
tains those terms proportional to (— /i"B + d"E), and 
T~Lfw2 contains those terms proportional to (^"E + (i"B). 
Focusing on the term proportional to (— yit'B + d'E), 
namely, [{J3^^)(n=2A),^i] in Eq. ^. For n = 2, we 
have 






(98) 
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where W is given in Eq. (|46)) and W^ is defined as 



-4c2^s • n(-^"B + d"E) • n, 



(99) 



which is proportional to electric and magnetic fields. By 

)(4),^b] 



using Eq. ([M)) . the term [{pQ.^)u\,Q.^] can be written as 



that the second term of Eq. (|100p is collected in h''^\ not 
/!(■*'. It can be shown that the only term that contributes 
toS-n(-/^"B + d"EVn|n|2 is [$0'^'^\0'-'^^] that is the 
correction term in h'^'^^ (see Eq. (jAlSp V Using O'^^ in 
Eq. ([94]), we have 
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It can be shown that [77, ai\ — —2/31]^, and thus we obtain 
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= P^^l + ^c4/3S • n(-M"B + d"E) • ninp. 



Consider those terms proportional to (— /^"B + (i"E): one 
comes form Eq. ([M|) in the corresponding term and the 
other is obtained form fi^;, Eq. dTOO]) and Eq. (fT02)) . Us- 
ing the definitions of ^ = Yi/mc (see Eq. ([TO))), /it" = -Eg/i' 
and d" = i^gd', one obtain 



^FWi = :t^W-^ + /3S ■ (-/i'B + d'E 
c'' /16 16 
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The first term of Eq. (|100p is just one of the 4th order 

r(4) 
'FW 



terms of -ffpw shown in Sec. |lTl It is important to note 



For (7e = 2 and ge — 2, T^fwi vanishes. We now trans- 
form ^ in Eq. pOSp in terms of the boost velocity (3. By 
using the transformation between ^ and /3 [see Eq. (|70|) ] , 
it can be shown that {\ - ||$|2) |^|2 « y|^|/3|2. Sub- 
stituting AMM coefficient ^' = (^ - l) ^ and AEM 
coefficient d 



' _ /9e 



2 1) 2!^ iiito Eq. (fT03)) . we find that 
Hpwi can be written as 
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On the other hand, the terms proportional to {fJ,"E + 
d"B) correspond to {floMl} = 2c^S -H x (/^"E + d"B). 
Similar to the derivation for "Hpwij we cohect ah terms 
proportional to {i7o, J7^} and obtain 



nFW2 = moM}^(i~lm' + lm' 



(105) 



Using /i" = Egfi' and d" ~ Egd' , we have 

Hfw2 =(1- ^m' + ^l^l^") S-^x (//E + d'B). (106) 
I 



By using Eq. ([TO)), we find that 

Hfw2 = S • /3 X (^'E + d'B) 

= (|-l)/x:^-(3xE) 



(107) 
1) < • X B). 



To focus on the interaction of dipole moments and exter- 
nal fields, we have to combine ilspin, "Hfwi and 'Hfw2 
together. After a straightforward calculations, we find 
that 



^spin + ^FWl + 'HfW2 
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Eq. (|108p is in agreement with Eq. (IM|) when the replace- 
ment 7^-7 and the duality transformation for elec- 
tromagnetic fields are used. Without magnetic charge, 
Eq. (|108p coincides with Eq. (^5)) for arbitrary values of 
Qe- The dual part of the TBMT equation for spin is also 
obtained. 

In short, the whole derivations in this section have as- 
sumed that electromagnetic fields are static and homo- 
geneous. Therefore, we show that up to terms of the 7th 
order in 1/ Eg, the FW transformation including anoma- 
lous dipole moments coincides with the TBMT equation 
for the spin-1/2 particle with arbitrary g^ and gg. 



VI. CONCLUSIONS AND DISCUSSION 

To investigate the low-energy limit of the relativistic 
quantum theory for a spin-1/2 charged particle, which 
is described by the Dirac equation, we perform a series 
of successive FW transformations on the Dirac Hamilto- 
nian up to terms of the 7th order in 1/Eg. Assuming the 
electromagnetic fields are static and homogeneous, and 
taking care of the relation between the kinematic mo- 
mentum n used in the Dirac Hamiltonian and the boost 



velocity (3 used in the TBMT equation, we show that 
the resulting FW transformation of the Dirac Hamilto- 
nian is in agreement with the classical orbital Hamilto- 
nian -fforbit plus the TBMT Hamiltonian ifspin with the 
gyromagnetic ratio ge being equal to 2. Through electro- 
magnetic duality, this can be generalized for a spin-1/2 
dyon, which has both electric and magnetic charges and 
thus possesses both intrinsic magnetic dipole moment fj,'^ 
and intrinsic electric dipole moment fi'p. 

To affirm the consistency between the low-energy limit 
of the relativistic quantum theory and the classical coun- 
terpart to a broader extent, we consider the relativistic 
quantum theory for a spin-1/2 dyon with arbitrary val- 
ues of the gyromagnetic and gyroelectric ratios, which is 
described by the Dirac-Pauli equation, namely, the Dirac 
equation with augmentation for AMM and AEM. Up the 
7th order in 1/Eg again, we show that the FW transfor- 
mation of the Dirac-Pauli Hamiltonian is also in accord 

with i^orbit + -f^spin- 

Many phenomena regarding spin dynamics have been 
observed and can be explained by the TBMT equation. 
These include the anomalous Zeeman effect, spin-orbit 
interaction, Thomas precession and chan ge r ate of the 
longitudinal polarization (see Sec. 11.11 of n^ for a brief 
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review). The TBMT equation is however derived merely 
by the requirement of covariance without invoking any 
first principles. By studying the FW transformation of 
the Dirac Hamiltonian and the Dirac-Pauli Hamiltonian, 
we have shown that the TBMT equation as a phenomeno- 
logical formula is in fact supported by the first principle 
of the fundamental relativistic quantum theory as a low- 
energy limit. (The relativistic quantum theory further 
requires the spin to be quantized as s = her/ 2; this re- 
sult cannot be obtained at the phenomenological level.) 
Therefore, the correspondence principle is again shown 
to be established. 

By far, the agreement between the Dirac/Dirac-Pauli 
equation and the orbital equation plus the TBMT equa- 
tion is only proven up to the 7th order in 1/ Eg. Further 
research is needed to investigate the FW transformation 
to higher orders and a generic expression for the FW 
transformation at an arbitrary order could be obtained 
by mathematical induction. If this is the case, performing 
successive FW transformations ad infinitum is expected 
to yield the result in precise agreement with the orbital 
equation plus the TBMT equation. 

Furthermore, the assumption of static and homoge- 
neous fields can be released. In time-varying and/or in- 
homogeneous fields, the TBMT equation has to be gener- 
alized to allow gradient force terms like (/i.^- V)B and the 
FW transformation of the Dirac/Dirac-Pauli Hamilto- 
nian shall yield the corresponding terms accordingly. The 
gradient force terms should not be missing, as (/x„j • V)B 
is used in the Stern-Gerlach experiment to separate spin- 
up and spin-down particles. Furthermore, the detailed 
investigation for the dipole moments interacting with the 
time- variation of electromagnetic fields may predict new 
physics. However, the calculation for the FW transfor- 
mation will be much more complicated if the fields are 
non-static and inhomogeneous. 
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Appendix A: Foldy-Wouthuysen transformation 



In this appendix, we expand the Dirac Hamiltonian up 
to terms of the 7th order in 1/Eg with Eg ~ 2mc^. The 
Dirac Hamiltonian can be separated in to two parts. One 
is the even operator denoted as $7^;, which commutes with 
/3, and the other is odd operator JIq, which anti-commutes 
with 0: 



{/3,fio} = 0. 
The Dirac Hamiltonian can be written as 



(Al) 



(A2) 



where ^Ie — ecf) + e(j) and flo — col ■ 11. The kinetic mo- 
mentum nisll = p— -A— -A. The first transformation 

^ c c 

operator can be written as 
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The Dirac Hamiltonian under the unitary transformation 
Ui can be written as 



J 



Hi = UiHU^^ ^^Eg + nE + Y, 



1 

Ef: 



(n + l) 

r 



^^)[mo)^n),no] + ^[m,)^n) 



,^e],, 



(A4) 



where the subscript n at (/3i7)(„) is defined as, for ex- as 
ample, [(;3f7„)(3), f]„] = [pQoAP^odP^c^E]]]- Equa- 
tion (jA4p can be again separated into odd and even parts. 
The even part of Eq. (|A4[) denoted as h can be written 



h = /i(") + 



n=l 9 



(A5) 
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where /i*^"^ = ^k and 



Therefore, Eq. (jA4p beconies 



^(«=1.3,5,-) _ 

y-nl (n + 1)! 



[(/3^o)(n),^o] 



(A6) 



E\ 



The odd part of Eq. (|A4p denoted as O can be written 
as 



Hi = '-Eg + h + O, 



(A9) 



'^^E 



©(^ 



EJ} 

ri=l 9 



(A7) 



where 



n! 

0(«=2,4,6.---) ^ ( ]_ _ 



(n+1)! 



[(/31]o)(„),n„]. 



(A8) 



where /i contains those terms with only even matrices and 
O contains only odd matrices. The second transforma- 
tion denoted as U2 — exp(S'2), where S2 is 6*2 — pO/Eg. 
We have 



H2 = U2HiU^' - !^Eg 



n=l 9 



;Jy-(;^)[(/30)(„),0] + ^[(^0)(„),/.] 



(3 



(AlO) 



E, + h' + O' 



where h' and C are the new even and odd parts, respectively, of the right hand side of the first equality. The even 
part of Eq. (JAIOI) denoted as h' can be written as 



?. ifs-hrhyi)!"^"'-"""! 



h' = h+ I 



n=l,3,5- 



(n+1)! 

n=2,4,6- a "■■ 



°° h'irn.) 
m=\ 3 



The odd term O' in Eq. (IA10[) is given by 



(All) 



O' = I 



y -i- 



1 



(n + 1)! 
1 1 



[{PO\n),0] 



n=l,3,5-- S 



\{iiO)(n)M 



E 



Ql(ra) 



m=3 9 



(A12) 



Firstly, we note that O'^^^ and C(^) in Eq. (jA12p vanish as we have m starting from 3. For the former result, the reason 
is that the lowest order of the unprimed odd term O is 1, and thus in the second line of the first equality in Eq. (jA12p . 
the primed odd term C is at least of the 2nd order. On the other hand, the explicit form of C(^) can be written as 
C"(2) = [/30(i\/i(")]. However, O^^) is given by Eq. (|M| for n = 1, namely, O^-^l = $[flo,flE] = icha ■ (eE + eE), 
which commutes with h^^' — fi^; — e(j) + ecj). If we perform the third transformation which is 53 = $0' /Eg at H2 — 
^ + h' + O', we will obtain a new even term h" as well as the new odd term O": H3 = 11:^11211:;^^ = ^Eg + h" + O" . 
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The even term h" is given by 



?. i(^-7;rh)i)i<^°''<», 



h" = h' +{ 



n=l,3,5' 



O'] 



n=2,4,6- ^S ^■ 



m=l a 



The odd term O" is given by 



o" = l 



E 



n=2,4,6- 



± 1 '— 

El' V"! (n + iy. 



[W0'\„^,0'] 



E i^;!iPo')H,^'] 



n=1.3.5...^."" 



E 



0"(™) 



(A13) 



(A14) 



The odd term O" starts form m = 4. Obviously, O'^^* is zero because C*-'"-' starts from m = 3. This can also be seen 
as follows. The explicit form of the next three terms of C)"(™) are 






(A15) 



Because O'^^^ and O'^^^ are zero, C'^^^ and ©"(^^ vanish. It can be shown that 0"'^^'> does not vanish. Using Eq. ((M3|) . 
/j"(™) for TO = 1, 2, ■ • ■ 6 are given by 



2!' 



(A16) 






(™\/l'(")]]. 



(f+m=4) 
4 



(^+m+n=3) 
3 2 



/i"(6) = ;,'(6) + ( 1 - i ) J2 [$0'^'\0'^"'^] + ^ E Et^^'^"' [/3C^'^"', /i'^"^] 



(^+m=5) 



2! 

{l+m+n=4:) 



Because O'^-^^ and C'^^) are zero, we have /i"^'^' = h''-^^ and /i''^*) = h'^^\ since the constraints £ + to == 4 and 
£ + m = 5 imply (£, to) = {(1, 3), (2, 2), (3, 1)} and {i,m) = {(1,4), (2, 3), (3, 2), (4, 1)}, respectively The commutator 
[^0/(f)^0/(m)j jj^ ;j//(5) j^^jj ;j»(6) vanishes. Consider the term [/30'(^), [/JO'^™), ft'w]] in /i"(5) and /i"^^) subject to 
the constraints £ + m + n = 3 and £ + m + n — 4, respectively. For n = 0, we have {£,m) = {(1,2), (2, 1)} and 
(£,to) = {(1,3), (2,2), (3, 1)} and thus [pO'^^\ [/30'("), /i'("="']] vanishes in h"^^^ and h'^^^K For n = 1 and n = 2, the 
term [/30'(^\ [/JO'^™), /i'(")]] still vanishes. 
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Therefore, up to terms of the 7th order in 7^, we obtain an important result 

/i"(") = /i'("\ n = l,2,---6, 
and /i'(") (i.e., Eq. (jAlip ) is given by 

/l'(2)=/i(2), 



{l+m=A) (^+m+n=3) 









(A17) 



(A18) 



Therefore, in order to obtain the FW transformation up to l/i?I, we need to know four odd terms O^^^ for n = 1, 2, 3, 4. 
On the other hand, if we perform the transformation again by using S'4 — J30" / Eg, since C)"(™) starts form tti = 4 
(i.e. 0"(™) = for m = 1,2,3), the transformation S'4 does not change Eq. (|XT7)) . The resulting C'"'(™) wih start 
from at least m — 5. To bring the odd term to the 7th order, we need S5 = $0"'/Eg and 5*6 = $0"" /Eg. However, 
the two transformations also do not change the validity of Eq. (|A17|) . Therefore the resulting Hamiltonian can be 
written as 



iJpw — UpwHUrfy^ 



P 



FW 
6 



^. + E4'^ + °(v4)- 



(A19) 



By using Eqs. (jA6p . (IA18|) and (|A17p . after straightforward calculations it can be shown that 






H 



(2) 
FW 

(3) 

FW - p3 

a 



1 /W 



E^g\2 



L^-^nt + H^-D)'}, 



(4) 



1 / 1 

a 



FW - ^ [^[[^o,w].^o] ~ ^[v,nl]] , 



(A20a) 
(A20b) 
(A20c) 
(A20d) 



HFW = 4^\l^A(^^oh,^^o] + l E W0^'\0("^^] + l ^ ^[/30W,[/30(™), /.(")]] 



£,m=l 



e,m=ln=0 
(«+m+n=3) 



<w = ^\^mo\^)^^o] + l E [/3o''\o^'"^] + J E Et/^^^'M/so'-^/^^"']] 



£6 ] 720 



(^+m=5) 



(£+m+n=4) 



(A20e) 

(A20f) 
(A20g) 
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where V = [no, Ob] and W = [D, fio]. 



Appendix B: Derivation of Equation (|57b|) 



The expHcit form of V,oWflo [Eq. ([57]) ] plays an im- 
portant role in obtaining the correct coefBcient of each 
term in H^^. By definition, Qo = cd • 11 and W — 
[[fio,rj£;],rio] = [[ca-U,V],ca-n], where V = e0 + e0. 
In the derivation for Eq. ([57)) . electromagnetic fields are 
assumed to be homogeneous and static. This means 
that the terms involving gradient of fields are neglected. 
Therefore, we have 



where Eq. (jB3p is used. The third term of Eq. (|B4p can 
be written as 



fih 

= epqr'^plli ( — ^rilBi + 11^11, 



ih 



(E.n.fi, - ^iXiqBt) + ep,,Ep|n|2n,. 



(B7) 



-^^oWVLo = -2c2;iep9^£:,n,n^nj(d.i]pdj), (bi) 



The fourth term of Eq. (IB4p becomes 



where £q = [eEq + e-Eg) and 

W = -2c2?iS • (5 X n) 
is used. Furthermore, it can be shown that 

ih 
[11^,11^] = —SijkBk, 

were Bk = eB^ + eBk- By using [a^, ffp] = 2ieip 
^pqr^ipm — 5qm^ri ~ SgiSmi, Eq. (jBip Can bc Written as 



(B2) 



(B3) 



and 



^^pqr^ijm^m^j^ "^ ''^ J 



ih 



(B8) 



— £pqr{2B£j2eIlr ~ BilliY^r + Bm^m)'^p, 

where eij^HiHj — ^Bm and eijm^rje = ^ir^mi ~ ^iAnr 
are used. We note that there is a field £q in the Eq. (|B4p . 
By neglecting the product of fields EiBj, Eq. (IB4p with 
substitutions of Eqs. (US]), (HU, (lB7| and dlH]) becomes 






~2T£q{2iaqaj\Y\.\ lij — 2iarCijTiqIir^j 

6pgr ZjplljllT-il^ -r l^pqr^.ijm^m^p^^i^^r^^j }• 



(B4) 



By using agcij 
becomes 



-ieqjfTii, the first term of Eq. (jB4 






-2c^ft£,(-2e,jfE^|n|^nj- + epqrJ:p\n\^Ur 



-inpw. 



(B9) 



2id„d,|n|2n,- = 2i|n|^n„ - 2e„,^s^|n|^nj. (B5) 



j;guj| 



tqjiZ^i] 



and we have 



On the other hand, the second term of Eq. (|B4p can be 
written as 



fioWfio 



-c^inPw. 



(BIO) 



2idrQ;wn„nrn, 



'2inq\n\ 



2ih 



^iUqBe, (B6) 
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